Functional significance may underlie the taxonomic utility of single amino acid substitutions in conserved proteins

We hypothesized that some amino acid substitutions in conserved proteins that are strongly fixed by critical functional roles would show lineage-specific distributions. As an example of an archetypal conserved eukaryotic protein we considered the active site of ß-tubulin. Our analysis identified one amino acid substitution—ß-tubulin F224—which was highly lineage specific. Investigation of ß-tubulin for other phylogenetically restricted amino acids identified several with apparent specificity for well-defined phylogenetic groups. Intriguingly, none showed specificity for “supergroups” other than the unikonts. To understand why, we analysed the ß-tubulin Neighbor-Net and demonstrated a fundamental division between core ß-tubulins (plant-like) and divergent ß-tubulins (animal and fungal). F224 was almost completely restricted to the core ß-tubulins, while divergent ß-tubulins possessed Y224. Thus, our specific example offers insight into the restrictions associated with the co-evolution of ß-tubulin during the radiation of eukaryotes, underlining a fundamental dichotomy between F-type, core ß-tubulins and Y-type, divergent ß-tubulins. More broadly our study provides proof of principle for the taxonomic utility of critical amino acids in the active sites of conserved proteins

Disease prevention versus data privacy : using landcover maps to inform spatial epidemic models

The availability of epidemiological data in the early stages of an outbreak of an infectious disease is vital for modelers to make accurate predictions regarding the likely spread of disease and preferred intervention strategies. However, in some countries, the necessary demographic data are only available at an aggregate scale. We investigated the ability of models of livestock infectious diseases to predict epidemic spread and obtain optimal control policies in the event of imperfect, aggregated data. Taking a geographic information approach, we used land cover data to predict UK farm locations and investigated the influence of using these synthetic location data sets upon epidemiological predictions in the event of an outbreak of foot-and-mouth disease. When broadly classified land cover data were used to create synthetic farm locations, model predictions deviated significantly from those simulated on true data. However, when more resolved subclass land use data were used, moderate to highly accurate predictions of epidemic size, duration and optimal vaccination and ring culling strategies were obtained. This suggests that a geographic information approach may be useful where individual farm-level data are not available, to allow predictive analyses to be carried out regarding the likely spread of disease. This method can also be used for contingency planning in collaboration with policy makers to determine preferred control strategies in the event of a future outbreak of infectious disease in livestock

Patterns of Genome Evolution among the Microsporidian Parasites Encephalitozoon cuniculi, Antonospora locustae and Enterocytozoon bieneusi

Microsporidia are intracellular parasites that are highly-derived relatives of fungi. They have compacted genomes and, despite a high rate of sequence evolution, distantly related species can share high levels of gene order conservation. To date, only two species have been analysed in detail, and data from one of these largely consists of short genomic fragments. It is therefore difficult to determine how conservation has been maintained through microsporidian evolution, and impossible to identify whether certain regions are more prone to genomic stasis.Here, we analyse three large fragments of the Enterocytozoon bieneusi genome (in total 429 kbp), a species of medical significance. A total of 296 ORFs were identified, annotated and their context compared with Encephalitozoon cuniculi and Antonospora locustae. Overall, a high degree of conservation was found between all three species, and interestingly the level of conservation was similar in all three pairwise comparisons, despite the fact that A. locustae is more distantly related to E. cuniculi and E. bieneusi than either are to each other.Any two genes that are found together in any pair of genomes are more likely to be conserved in the third genome as well, suggesting that a core of genes tends to be conserved across the entire group. The mechanisms of rearrangments identified among microsporidian genomes were consistent with a very slow evolution of their architecture, as opposed to the very rapid sequence evolution reported for these parasites

Gene transfer: anything goes in plant mitochondria

Parasitic plants and their hosts have proven remarkably adept at exchanging fragments of mitochondrial DNA. Two recent studies provide important mechanistic insights into the pattern, process and consequences of horizontal gene transfer, demonstrating that genes can be transferred in large chunks and that gene conversion between foreign and native genes leads to intragenic mosaicism. A model involving duplicative horizontal gene transfer and differential gene conversion is proposed as a hitherto unrecognized source of genetic diversity

A broad distribution of the alternative oxidase in microsporidian parasites

Microsporidia are a group of obligate intracellular parasitic eukaryotes that were considered to be amitochondriate until the recent discovery of highly reduced mitochondrial organelles called mitosomes. Analysis of the complete genome of Encephalitozoon cuniculi revealed a highly reduced set of proteins in the organelle, mostly related to the assembly of ironsulphur clusters. Oxidative phosphorylation and the Krebs cycle proteins were absent, in keeping with the notion that the microsporidia and their mitosomes are anaerobic, as is the case for other mitosome bearing eukaryotes, such as Giardia. Here we provide evidence opening the possibility that mitosomes in a number of microsporidian lineages are not completely anaerobic. Specifically, we have identified and characterized a gene encoding the alternative oxidase (AOX), a typically mitochondrial terminal oxidase in eukaryotes, in the genomes of several distantly related microsporidian species, even though this gene is absent from the complete genome of E. cuniculi. In order to confirm that these genes encode functional proteins, AOX genes from both A. locustae and T. hominis were over-expressed in E. coli and AOX activity measured spectrophotometrically using ubiquinol-1 (UQ-1) as substrate. Both A. locustae and T. hominis AOX proteins reduced UQ-1 in a cyanide and antimycin-resistant manner that was sensitive to ascofuranone, a potent inhibitor of the trypanosomal AOX. The physiological role of AOX microsporidia may be to reoxidise reducing equivalents produced by glycolysis, in a manner comparable to that observed in trypanosome

Unexpected mitochondrial genome diversity revealed by targeted single-cell genomics of heterotrophic flagellated protists

This is the author accepted manuscript. The final version is available from Nature Research via the DOI in this recordData availability: Complete mtDNA sequences assembled from this study are available at GenBank under the accession numbers MK188935 to MK188947, MN082144 and MN082145. Sequencing data are available under NCBI BioProject PRJNA379597. Reads have been deposited at NCBI Sequence Read Archive with accession number SRP102236. Partial mtDNA contigs and other important contigs mentioned in the text are available from Figshare at https://doi.org/10.6084/m9.figshare.7314728. Nuclear SAG assemblies are available from Figshare at https://doi.org/10.6084/m9.figshare.7352966. A protocol is available from protocols.io at: https://doi.org/10.17504/protocols.io.ywpfxdn.Code availability: The bioinformatic workflow is available at https://doi.org/10.5281/zenodo.192677; additional statistical analysis code is available at https://doi.org/10.6084/m9.figshare.9884309.Most eukaryotic microbial diversity is uncultivated, under-studied and lacks nuclear genome data. Mitochondrial genome sampling is more comprehensive, but many phylogenetically important groups remain unsampled. Here, using a single-cell sorting approach combining tubulin-specific labelling with photopigment exclusion, we sorted flagellated heterotrophic unicellular eukaryotes from Pacific Ocean samples. We recovered 206 single amplified genomes, predominantly from underrepresented branches on the tree of life. Seventy single amplified genomes contained unique mitochondrial contigs, including 21 complete or near-complete mitochondrial genomes from formerly under-sampled phylogenetic branches, including telonemids, katablepharids, cercozoans and marine stramenopiles, effectively doubling the number of available samples of heterotrophic flagellate mitochondrial genomes. Collectively, these data identify a dynamic history of mitochondrial genome evolution including intron gain and loss, extensive patterns of genetic code variation and complex patterns of gene loss. Surprisingly, we found that stramenopile mitochondrial content is highly plastic, resembling patterns of variation previously observed only in plants.Gordon and Betty Moore FoundationLeverhulme TrustDavid and Lucile Packard FoundationRoyal SocietyEuropean Molecular Biology OrganizationCONICYT FONDECYTGenome Canad

Extensive molecular tinkering in the evolution of the membrane attachment mode of the Rheb GTPase

Rheb is a conserved and widespread Ras-like GTPase involved in cell growth regulation mediated by the (m)TORC1 kinase complex and implicated in tumourigenesis in humans. Rheb function depends on its association with membranes via prenylated C-terminus, a mechanism shared with many other eukaryotic GTPases. Strikingly, our analysis of a phylogenetically rich sample of Rheb sequences revealed that in multiple lineages this canonical and ancestral membrane attachment mode has been variously altered. The modifications include: (1) accretion to the N-terminus of two different phosphatidylinositol 3-phosphate-binding domains, PX in Cryptista (the fusion being the first proposed synapomorphy of this clade), and FYVE in Euglenozoa and the related undescribed flagellate SRT308; (2) acquisition of lipidic modifications of the N-terminal region, namely myristoylation and/or S-palmitoylation in seven different protist lineages; (3) acquisition of S-palmitoylation in the hypervariable C-terminal region of Rheb in apusomonads, convergently to some other Ras family proteins; (4) replacement of the C-terminal prenylation motif with four transmembrane segments in a novel Rheb paralog in the SAR clade; (5) loss of an evident C-terminal membrane attachment mechanism in Tremellomycetes and some Rheb paralogs of Euglenozoa. Rheb evolution is thus surprisingly dynamic and presents a spectacular example of molecular tinkering

Stops making sense: translational trade-offs and stop codon reassignment

Background Efficient gene expression involves a trade-off between (i) premature termination of protein synthesis; and (ii) readthrough, where the ribosome fails to dissociate at the terminal stop. Sense codons that are similar in sequence to stop codons are more susceptible to nonsense mutation, and are also likely to be more susceptible to transcriptional or translational errors causing premature termination. We therefore expect this trade-off to be influenced by the number of stop codons in the genetic code. Although genetic codes are highly constrained, stop codon number appears to be their most volatile feature. Results In the human genome, codons readily mutable to stops are underrepresented in coding sequences. We construct a simple mathematical model based on the relative likelihoods of premature termination and readthrough. When readthrough occurs, the resultant protein has a tail of amino acid residues incorrectly added to the C-terminus. Our results depend strongly on the number of stop codons in the genetic code. When the code has more stop codons, premature termination is relatively more likely, particularly for longer genes. When the code has fewer stop codons, the length of the tail added by readthrough will, on average, be longer, and thus more deleterious. Comparative analysis of taxa with a range of stop codon numbers suggests that genomes whose code includes more stop codons have shorter coding sequences. Conclusions We suggest that the differing trade-offs presented by alternative genetic codes may result in differences in genome structure. More speculatively, multiple stop codons may mitigate readthrough, counteracting the disadvantage of a higher rate of nonsense mutation. This could help explain the puzzling overrepresentation of stop codons in the canonical genetic code and most variants